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Abstract 

We consider one-dimensional Schrodinger-type operators in a bounded 
interval with non-self-adjoint Robin-type boundary conditions. It is well 
known that such operators are generically conjugate to normal operators 
via a similarity transformation. Motivated by recent interests in quasi- 
Hermitian Hamiltonians in quantum mechanics, we study properties of 
the transformations in detail. We show that they can be expressed as 
the sum of the identity and an integral Hilbert-Schmidt operator. In the 
case of parity and time reversal boundary conditions, we establish closed 
integral-type formulae for the similarity transformations, derive the simi- 
lar self-adjoint operator and also find the associated "charge conjugation" 
operator, which plays the role of fundamental symmetry in a Krein-space 
reformulation of the problem. 
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1 Introduction 



Let us consider the m-sectorial realization H of the second derivative operator 

ip i— > -ip" (1.1) 

in the Hilbert space % := L 2 (—a, a), with a > 0, subjected to separated, Robin- 
type boundary conditions 

i>'{±a) +c±4>{±a) = (1.2) 

where c± are arbitrary complex numbers. The operator H is self-adjoint if, 
and only if, the constants c± are real. The present paper is concerned with 
the existence and properties of similarity transformations of H to a normal or 
self-adjoint operator in the non-trivial case of non-real c±. 

The similarity to the normal (respectively, self-adjoint) operator is under- 
stood as the existence of a bounded operator Q with bounded inverse such that 

h:= Of/XT 1 (1.3) 

is normal (respectively, self- adjoint). We remark that this concept is equivalent 
to the existence of a topologically equivalent inner product in H with respect 
to which H is normal (respectively, self-adjoint). In addition to results on the 
general structure of the similarity transformations, modified inner products, and 
transformed operators, we present explicit closed formulae for these objects in 
special cases of boundary conditions. 

The operators of the type (|l.ip - (|1.2p have been studied from many aspects 
and there exist a large number of known results; we particularly mention the 
classical monograph of Dunford and Schwartz [TTJ Chapter XIX. 3]. Recent years 
brought new motivations and focused attention to some aspects of the problem 
which attracted little attention earlier. 

As an example, let us mention that one-dimensional Schrodinger operators 
with non-Hermitian boundary conditions of the type (jl.2p were used as a model 
in semiconductor physics by Kaiser, Neidhardt and Rehberg [18j . In their paper 
the imaginary parts of the constants c± are required to have opposite signs such 
that the system is dissipative. The authors find the characteristic function 
of the operators, construct its minimal self-adjoint dilation and develop the 
generalized eigenfunction expansion for the dilation. See also [161 117] for further 
generalizations. Here the main idea of using non-self-adjointness comes from 
embedding a quantum-mechanically described structure into a macroscopic flow 
and regarding the system as an open one. 

However, the principal motivation of the present work is the possibility of 
giving a direct quantum-mechanical interpretation of non-Hermitian operators 
which are similar to self-adjoint ones [29) ■ The most recent strong impetus to 
this point of view comes from the so-called PT-symmetric quantum mechanics. 
Here the reality of the spectrum of a class of non-Hermitian operators - caused 
by certain symmetries rather than self-adjointness - suggests their potential 
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relevance as quantum-mechanical Hamiltonians; see the review articles [H |2"T]. 
It has been confirmed during the last years that it is indeed the case provided 
that the similarity transformation to a self-adjoint operator can be ensured. 
However, it is a difficult task. 

Motivated by the lack of rigorous results, the authors of [2T] introduced 
a simple non-Hermitian 'PT-symmetric operator of the type (ll.ip - (|1.2l) and 
wrote down a closed formula for the (square of the) similarity transformation 
(see also [501 [35]). Let us also mention that the importance of (not only) PT- 
symmetric version of (|1.1I) - (|1.2I) in quantum mechanical scattering has been 
recently established in [T3] . 

The present paper can be regarded as a step further. In addition to consider- 
ing more general situations of larger classes of boundary conditions and similar- 
ity to normal operators, we provide an alternative and more elegant (integral- 
kernel) formulae for the similarity transformations in the PT-symmetric sit- 
uation. Moreover, we also give a remarkably simple formula for the similar 
self-adjoint operator in this case. Finally, we succeed in finding the so-called 
C-operator in a closed form, which plays the role of fundamental symmetry in 
a Krein-space reformulation of the problem. 

The paper is organised as follows. In Section [5J we give a precise definition 
of the operator H, summarize its known properties and recall the general con- 
cepts of quasi- Hermitian, 'PT-symmetric, and C-symmetric operators. Our main 
results about the universal structure of the similarity transformations can be 
found in Section [3] In Section [4] we show how these can be applied to particular 
(PT-symmetric) classes of boundary conditions and we present some explicit 
constructions of the studied objects. In Section[5]we discuss how the results can 
be extended to bounded and even second-order perturbations of H . Our final 
Section [5] presents a series of concluding remarks. 



2 Preliminaries 

We start with recalling general properties of H and concepts of similarity trans- 
formations in Hilbert spaces. 

2.1 Definition of the operator H 

The standard norm in our Hilbert space % = L 2 (—a,a) is denoted by || ■ ||. 
The corresponding inner product is denoted by (•,•) and it is assumed to be 
antilinear in the first component. 

We consider the m-sectorial realization H of the operator (|1.1[) subjected 
to the boundary conditions (|1.2I) as the operator associated on H with the 
quadratic form 

tHM:=ll^l| 2 + c + |^(a)| 2 - C -l^(-«)| 2 , 
Dom(iff) := W h2 (-a,a). 
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Note that the boundary terms are well defined because of the embedding of 
the Sobolev space IV 1 ' 2 (— a, a) in the space of uniformly continuous functions 
C°[—a,a], An elementary idea of the proof of the embedding can be also used 
to show that the boundary terms represent a relatively bounded perturbation 
of the form associated with the Neumann Laplacian (i.e., c± = 0). Since the 
Neumann form is clearly non-negative and closed by definition of the Sobolev 
space, we know that tn is a closed sectorial form by a standard perturbative 
argument QH Sec. VI. 1.6]. 

By the representation theorem [191 Thm. VI. 2.1] and an elementary version 
of standard elliptic regularity theory, it is easy to see that 

Htb = -ib", 

(2 2) 

Dom (H) = {ip g TV 2 ' 2 (-a, a) : ik'(±a) + c±ip(±a) = 0}. 

We refer to [TTJl Ex. VI. 2. 16] for more details. The operator definition (|2.2j) 
gives a precise meaning to (|1.1|) - (|1.2|) . 



2.2 Dirichlet and Neumann boundary conditions 

This subsection is mainly intended to collect some notation we shall use later. 

We have already mentioned that the special choice c± — gives rise to 
the Neumann Laplacian — A at on T~L. The Dirichlet Laplacian — Ad on H can 
be considered as the other extreme case by formally putting c± — +oo. It 
is properly defined as the second derivative operator (jl.ip with the operator 
domain Dom (~A D ) := W 2 - 2 (-a 1 a) n W^ 2 (-a, a). 

The spectrum of the Dirichlet and Neumann Laplacians in our one-dimen- 
sional situation is well known: 

a(-A D ) = {k 2 n }%L l , 

„ with k n := — — . 

*(-A N ) = {k 2 n }%L , 2a 

The corresponding eigenfunctions are respectively given by 

cos k n ( x + a) if n > 1 . 

(2.3) 

To simplify some expressions in the sequel, we extend the notation by \o '■= 0- 
Next we introduce a "momentum" operator p and its adjoint p*: 



n 1 n I if n = , 

Xn( x ) '■= -i=smk n (x + a), Xn( x ) 
yja 



pip := —itp', p*ib — —iib', 

Dom(p) := Wo' 2 (-a,a), Dom(p*) = W 1>2 (-a,o). 



The following identities hold: 



D jt, N ■ * N _ i D 

l PXn K n Xn ) Xn — ^nXn > 

-A D = p*p, - Ajv = pp* ■ 



(2.4) 



(2.5) 
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The resolvents (—Ad — k 2 ) 1 , (— An — k 2 ) 1 act as integral operators with 
simple kernels (Green's functions) Q 1 ^ and Q^, respectively: 

k _ - sin(k(x + a)) sm(k(y-a)) 

Qd{X ' V) - ' * <V ' (26) 

g^^y)= " c ° s(fc ^ +a) ^°y ( "" a)) ' 

fc sm(ifca) 

with x, y exchanged for x > y. Here the spectral parameter k 2 is supposed to 
belong to the resolvent set of the respective operator. 
For k = 0, the kernel of (— Ad)" 1 simplifies to 

no t \ {x + a)(a-y) (t) 7 , 

GD{x,y) = Y a > x<y, (2.7) 



with x, y exchanged for x > y. The resolvent of — A^v does not exist for k = 0, 
of course, but one can still introduce the reduced resolvent (— A^) of the 
Neumann Laplacian with respect to the eigenvalue (see [T9l Sec. III. 6. 5] for the 
concept of reduced resolvent). From the point of view of the spectral theorem: 



(-A^-^g-^x^M- (2.8) 

n=l " 

The corresponding integral kernel Gjj{x,y) can be obtained by taking the limit 
k — > of the regularized expression G| r (x, y) + k~ 2 Xo ( x )Xo (v)- We find 

^) = ^ + ^-f> *<V> ^ 

with x, y exchanged for x > y. 
Finally, we introduce operators 

oo 

J'~J2 C nXn(X^-), oe{D,N}, (2.10) 

where C n are positive numbers satisfying 

< mi < C n < m 2 < oo ( 2 -ll) 

for all n > 0, with given positive m 1 ,m2. The sum in the definition (|2.10l) . as 
well as all other analogous expressions in the following, are understood as limits 
in the strong sense. 

2.3 General properties of H 

Now we are in a position to recall some general properties of the operator H . 
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Proposition 2.1 (General known facts). 

(i) H is m-sectorial. The adjoint operator H* is obtained by taking the com- 
plex conjugation of c± in the boundary conditions (|1.2[) . 

(ii) H forms a holomorphic family of operators of type (B) with respect to the 
boundary parameters c± . 

(iii) The resolvent of H is a compact operator. 

(iv) H is a discrete spectral operator. 

(v) If all eigenvalues are simple, then H is similar to a normal operator. If 
the spectrum of H is in addition real, then H is similar to a self-adjoint 
operator. 

We have already shown that H is m-sectorial as the operator associated with 
the closed sectorial form (I2.2[) . The rest of the claim (i) follows by the fact that 
the adjoint operator is associated with the adjoint form t* H ((j),ip) := tij(ip,(j>) 
(cf. OH Thm. VI.2.5]). Property (ii) follows from §X?2§ as well if we recall that 
the boundary terms represent a relatively bounded perturbation of the form 
associated with the Neumann Laplacian and the relative bound can be made 
arbitrarily small {cf. [HI Sec. VII.4.3]). This also proves (iii) as a consequence 
of the perturbation result [HJJ Sec. VI. 3. 4]. The proof of (iv) is contained in [TT1 
Chapter XIX. 3]. Property (v) is a consequence of (iv). 

The similarity to a normal operator can be equivalently stated as the Riesz 
basicity of the eigenvectors of H . This property is shared by all second derivative 
operators with strongly regular boundary conditions, see [26] . Using the notion 
of spectral operator, this has been investigated in [IT] as well. 

Although the eigenvalues of H are generically simple, degeneracies may ap- 
pear. However, the only possibility are the eigenvalues of algebraic multiplicity 
two and geometric multiplicity one. In this case, operator H cannot be sim- 
ilar to a normal one, nevertheless, the eigenvectors together with generalized 
eigenvectors still form a Riesz basis. 

Now we turn to symmetry properties of H. 

Definition 2.1 (T'T-symmetry) . We say that H is VT -symmetric if 

[PT,H}=0, (2.12) 

where 

(^)W:=^N, (TV)(z) (2-13) 

It should be stressed that VT is an antilinear operator. The commutator 
relation (I2.12[) means precisely that (VT)H C H(V7~), as usual for the commu- 
tativity of an unbounded operator with a bounded one (cf. [HI Sec. III. 5. 6]). In 
the quantum-mechanical context, V corresponds to the parity inversion (space 
reflection), while T is the time reversal operator. 
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Definition 2.2 (S'-self-adjointness). We say that H is S '-self -adjoint if the re- 
lation H = S~ l H* S holds with a boundedly invertible operator S. 

We will use this concept in a wide sense, with S being either linear or anti- 
linear operator. If S is a conjugation operator (i.e. antilinear involution), then 
our definition coincides with the concept of J-self-adjointness [T21 Sec. III. 5]. 

While Definition ^. 2l is quite general, Definition ^. ll makes sense for operators 
in a complex functional Hilbert space only. In our case, we have: 

Proposition 2.2 (Symmetry properties). 

(i) H is T -self-adjoint. 

(ii) H is V -self-adjoint if and only if, c_ = — c+. 

(iii) H is VT - symmetric if, and only if, c_ = — c+ . 

Property (ii) coincides with the notion of self-adjointness in the Krein space 
equipped with the indefinite inner product (-,V-). It is also referred to as V- 
pseudo-Hermiticity in physical literature (see, e.g., [27]). 

It follows from Proposition 12. 21 (i) that the residual spectrum of H is empty 
(cf. Corol. 2.1]). Alternatively, it is a consequence of Proposition 12.11 fiii) . 
which in addition implies that the spectrum of H is purely discrete. 

We denote the (countable) set of eigenvalues of H by {A„}^L and the corre- 
sponding set of eigenfunctions by {V'nj^o- Similarly, let {A„}^ and {0 n }^o 
be the set of eigenvalues and eigenfunctions of the adjoint operator H* . That is 

Htp n = X„ip n , H*4> n = X l <t>n- (2.14) 

Eigenfunctions tp„ and 4>m corresponding to different eigenvalues, i.e. X n ^ \ rn , 
are clearly orthogonal. Solving the eigenvalue equation for H in terms of sine 
and cosine functions, it is straightforward to reduce the boundary value problem 
to an algebraic one. 

Proposition 2.3 (Spectrum). The eigenvalues X n = i„ of H are solutions of 
the implicit equation 

sin(2a/)(c_c+ + l 2 ) + (c_ - c+)l cos(2aZ) = . (2.15) 

The corresponding eigenfunctions of H and H* respectively read 



1 / c _ 
■)p n (x) = A n —= cos(l n (x + a)) — — sm(l n (x + a)) 
V a V '» 

If— — — 

<j) n {x) = cos(Z n (cc + a)) — = sm(l n (x + a)) 

V a V In 



(2.16) 



// all eigenvalues are simple, ip n can be normalized through the coefficients A n 
in such a way that (i^ n ,4'm) = S nm . 
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The spectrum of H has been described more explicitly for the PT-symmetric 
case. First of all, as a consequence of the symmetry, we know that the spectrum 
is symmetric with respect to the real axis. In the following proposition we 
summarize more precise results obtained in |21[ 122] . 

Proposition 2.4 (PT-symmetric spectrum). Let c± —ia± f3, with a, j3 € K. 
1. If j} = then all eigenvalues of H are real, 

A = a 2 , A„ = kl, ne N. (2.17) 
The corresponding eigen] 'unctions of H and H* respectively read 



Mx) = 4= eia[X+a) ' M*) = *Z W + if X» (*)• 
V2a k n 



(2.18) 



If a 7^ fe„ /or every n 6 N, i/ien a// t/ie eigenvalues are simple and choosing 



fyp 2{aa \P2n k 2 
Ao:= ae ,/f , (2.19) 
sm(2aa) fc„ — a z 

we /iawe i/ie biorthonormal relations (ipmfim) = S nm - 

2. If {3 > 0, t/ien t/ie eigenvalues of H are real and simple. 

3. If f3 < 0, i/ien a/Z the eigenvalues are either real or there is one pair of 
complex conjugated eigenvalues with real part located in the neighborhood 
of a 2 + f3 2 . ' 

In any case, the eigenvalue equation (|2 . 15[) can be rewritten as 

{I 2 - a 2 - (3 2 ) sin(2aZ) - 2(31 cos(2ai) = 0. (2.20) 



2.4 Concept of the metric operator 

We recall the concept of metric operator (or quasi-Hcrmitian operators intro- 
duced in [9]), widely used in 'PT-symmetric literature. 

Definition 2.3 (Metric operator and quasi-Hermiticity). Bounded positive^ 
operator with bounded inverse is called a metric operator for H , if H is 0- 
self- adjoint. H is then called quasi- Hermitian. 

It is obvious that the quasi-Hermitian operator H is self-adjoint with respect 
to the modified inner product (-,-)e : — (")©')■ It is also not difficult to show 
that the metric operator exists if, and only if, H is similar to a self-adjoint 

1 A is positive if {/, Af) > for all / € H, f / 0. 
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operator. Moreover, since H has purely discrete spectrum, the metric operator 
can be obtained as 

oo 

e = ^c r 2 >„ (<£„,.), (2.21) 

n=0 

where cf> n are eigenfunctions of H* and C n are real constants satisfying (|2.11l) . 

The expression (|2.2ip illustrates a non-uniqueness of the metric operator 
caused by the arbitrariness of C n . The latter can be actually viewed as a 
modification of the normalization of functions 4> n . Choosing different sequences 
{C n }^ = o, we obtain all metric operators for H, cf. [3Tj 133] . 

It is important to stress that if we define an operator O by (|2.21[) . we find 
that such O is bounded, positive, and with bounded inverse whenever {<^n}^=o 
is a Riesz basis. Thus, by virtue of Proposition |2.11 (v). such a O exists if, 
and only if, all eigenvalues of H are simple. However, the 0-self-adjointness 
of H is satisfied if, and only if, the spectrum of H is real. Otherwise, only 
QHQ^H* = H*QH&- 1 holds, cf. [33], which is equivalent to the fact that H 
is similar to a normal operator. 

In the following, the operator is always defined by (12.21)) regardless if it 
is a metric operator for H in view of Definition 12.31 

It should be also noted that 0, as a positive operator, can be always decom- 
posed to 

= n*Q. (2.22) 

One example of such Q is obviously We shall take the advantage of some 
different decompositions of the type (|2.22[) later. It follows easily from Defini- 
tion 12.31 that the similar operator h defined by (|1.3j) with il given by (|2.22j) is 
self-adjoint if is a metric operator for H. If all eigenvalues of H are simple 
but no longer entirely real, h is (only) a normal operator. 

2.5 Concept of the C operator 

For PT-symmetric operators, the notion of C operator was introduced in [B] and 
formalized in [2] . It was observed in [25] and in many works after that paper that 
Krein spaces provide suitable framework for studying PT-symmetric operators. 
Indeed, PT-symmetric operators which are at the same time P-self-adjoint are 
in fact self-adjoint in the Krein space equipped with the indefinite inner product 
(•,P-). Recall that our operator H is P-self-adjoint if, and only if, it is PT- 
symmetric {cf. Proposition ^. 2p . 

Definition 2.4 (C operator). Assume that H is V -self-adjoint (cf. Proposi- 
tion \2.2\ ). We say that H possesses the property of ' C- symmetry, if there exists 
a bounded linear operator C such that [H,C] — 0, C 2 = I, and PC is a metric 
operator for H . 

Thus, from the point of view of metric operators, we can find the C operator 
as C := PO for satisfying (PO) 2 = /. Hence C-symmetry allows us to 
naturally choose a metric operator. Besides a possible physical interpretation 
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of C discussed in [51 3] , it appears naturally in the Krein spaces framework as 
pointed out in [23J [24] as a fundamental symmetry of the Krein space (H, (•, V-)) 
with an underlying Hilbert space (H 7 (-,VC-)), 

3 General results 

In this section we provide general properties of the metric operator O defined 
in (|2.2ip and its decompositions fl from (|2.22[) . 

Let {^ n }%Lo and {0„}^L O denote the set of eigenvectors of H and H* , 
respectively. We assume that tp n and (f) n form Riesz bases and that they are 
normalized in such a way that (xp n , <f> m ) — S rnn . In view of Propositions l2.11 12.31 
we know that this is satisfied if all the eigenvalues of H are simple, which is a 
generic situation. 

Let {en} $^=0 be any orthonormal basis of H. If all eigenvalues of H are 
simple, we introduce an operator O by 

oo 

Q:=J2 e n{<t>n,-). (3.1) 
n=0 

Clearly, Cl : ip n i-> e„. 

ft is defined only if all eigenvalues are simple, however, sometimes it is pos- 
sible to extend it by continuity, see examples in Section 01 Nonetheless, such 
Q is typically not invertible and the dimension of the kernel corresponds to the 
size of Jordan blocks appearing in the spectrum of H . 

Basic properties of fi are summarized in the following. 

Lemma 3.1. Let all eigenvalues of H be simple. Then Q is a bounded operator 
with bounded inverse given by 

oo 

n_1 = 2^n(en,->, (3-2) 

?i=0 

i.e. fl^ 1 : e n M> ip n . The adjoint of Q reads 

oo 

fl* = J2^(e n ,-)- (3-3) 

i.e. Q* : e„ t-¥ 4> n and = 0, where G is defined in (|2 . 2 1 1) with C n = 1. 

Furthermore, we show how the operator Q can be realized. 

Theorem 3.2. Let all eigenvalues of H be simple, can be expressed as 

fl = U + L, (3.4) 

where U :— X)n°=o e n(Xn ■> ')> U '■ Xn e «> * s a unitary operator, and L is a 
Hilbert- Schmidt operator. 
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Proof. At first we remark that it suffices to prove that ft = I + L for e n := 
Xn j where £ is Hilbcrt-Schmidt. More precisely, if we compose U from the 
claim and 7 + L, we obtain £1 in (13.41) since L = UL is Hilbert-Schmidt too. 
Thus, we consider this choice of e„ in the following. Furthermore, we put a := 
7r/2 for simplification of the formulae. This specific choice is in fact harmless, 
since we can easily transfer the results for different a using the isometry V : 
J L 2 (-vr/2,7r/2) -> L 2 (-a,a) defined by ip{x) i-> y|T^(|2i). 

The asymptotic analysis of eigenvalues of H in [TTJ proof of Lem. XIX. 3. 10] 
shows that 

l n = n+ C+ ~ C - +0(n~ 2 ), 

J" 1 _ v (3-5) 
A n ^ = fe* + Ac + C - J +0(r t - 1 ), 

and |Im(Z„)| is uniformly bounded in n. These formulae are valid except for a 
finite number 7V of eigenvalues. 

We set e„ := i n — k n = l n — n. Using elementary trigonometric identities, 
we rewrite the eigenfunctions <p n as follows 



<t>n(x) = Xn ( x ) cos ( £ n( x + a)) ~ Xn ( X ) sin ( £ n( x + a )) 

- = [Xn( x ) cos(e^(x + a)) + Xn ( x ) sin i^( x + a ))] ■ 

We further rewrite the cosine and sine functions in this expression as 

, , xx , 2 cos fe(x + a)) — 1 , 

cos (e n (a; + a)) = 1 + e„ =5 =: 1 + e n c n (x), 

. ,—, , ss _ sin (e^(x + a)) _ 
sin (e„(a; + a)) = e„ — =: e„ s n (x). 



(3.6) 



(3.7) 



Note that ||c„|| and ||s„|| are uniformly bounded in n because of the properties 
of e n . The building block Xn {<Pn, •) of ft then becomes 

X^{K, •> = X^ (x%, •> + <&tf<X? c„, •) - e n xZ(x5s n , ■> _^ § 

i (Xn (Xn ' ') + £ nXn (Xn C m ') + £ "Xn (X n s nj ')) • 

Taking the sum of Xn (^m ') as m <|3-1[) , we obviously get f2 = 7 + L. 

It remains to show that the Hilbert-Schmidt norm ||7||hs of L is finite. We 
will understand L as a sum L — Ljy + 7oo > where 

AT -1 oo 

L No ■■= Xnikr), ioo := Xnikr), (3.9) 

n=0 n=N 

and <^„ := </>„ — Xn ■ ^n is a finite rank operator, hence it is automatically 
Hilbert-Schmidt and it suffices to consider 7^ in the rest of the proof. We 
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estimate explicitly only one term in the expression for ||Xoo|IhS! * ne res ^ follows 
in a similar way: 



oo / oo 



y , ( y , £ nXn {Xn s niXp ) I ^ ] e mXm {x m s m, Xp ) ) 
p=0 \n=N m=N I 

oo oo -.00 

= E E K?\{x^s n ^)\ 2 <- n E M 2 KII 2 <oc. 



(3.10) 



a 

p=0 n=N n=N 



Here the first inequality follows by the Bessel inequality (after interchanging the 
order of summation, which is justified) and by estimating \n by its supremum 
norm. The asymptotic behavior of e n and the uniform boundedness of ||s n || are 
used in the last step. □ 

Corollary 3.3. Let all eigenvalues of H be simple. Then 

6:=n*n = I + K (3.11) 

coincides with defined in (|2.21l) with C n = 1. Here K is a Hilbert- Schmidt 
operator that can be realized as an integral operator with a kernel belonging to 
L 2 ((—a,a) x (—a, a)). 



Proof. The claim follows from Theorem 13.21 and the well-known facts that 
Hilbert-Schmidt operators are *-both-sided ideal in the space of bounded oper- 
ators and can be realized as integral ones, see [28j Thm.VI.23]. □ 

Remark 3.1. Slight modification of the definition of f2 and the proof of The- 
orem 13.21 yields the analogous result for operators O defined in (|2.21[) with 
arbitrary C n . It suffices to consider /„ :— C n e n instead of e n . The resulting 
form is 

e = J N + K, (3.12) 

where J N is defined in (|2 . 10[) and K is again a Hilbert-Schmidt operator. J N it- 
self, however, can be a sum of a bounded and a Hilbert-Schmidt operator, as 
we shall see in examples. 

Proposition 3.4. Let S be an open connected set in C 2 such that for all 
(c_,c + ) € S all eigenvalues of H are simple. Then and thereby are bounded 
holomorphic families in S with respect to parameters c± . 

Proof. We verify the criterion stated in [191 Sec. VII. 1.1]. We have proved 
already that £1 is bounded. It remains to show that (/, fig) is holomorphic for 
every /, g from a fundamental set of % that we choose as the orthonormal basis 
{ e n}^Lo- ( e mj^ e n) = (0mj e n) is holomorphic because 4> m is an eigenfunction 
of the operator H* , which can be viewed as a holomorphic family of operators 
of type (B) with respect to the parameters c±. □ 

Corollary 3.5. Assume the hypothesis of Proposition \3^4\ Then h :— flHfl^ 1 
is a holomorphic family of operators in S with respect to parameters c± . 
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Since the operator H is a holomorphic family of type (B), i.e. it is naturally 
defined via quadratic forms with the domain W 1,2 (— a, a) independent of the 
parameters c±, h is expected to possess a similar property. To prove it, we 
have to particularly show that the associated quadratic forms corresponding 
to different values of c± have the same domain, which is not guaranteed by 
Corollarv l3.5l To this end we analyse the quadratic form associated to h, where 
we set e„ := Xn m the definition of fi. 

Theorem 3.6. Let all eigenvalues of H be simple and let e n := \n * n CO}- 

Then = I + L and fl^ 1 = I + M , where L, M are Hilbert- Schmidt opera- 
tors, fi, 0,*, f2 _1 , (fi -1 )* are bounded operators on W 1 ' 2 (—a, a) andW 2 ' 2 (—a,a). 
Furthermore, the following estimates hold for all <fi £ W ' {— a, a) and arbitrary 

||(i»'|| 2 <C( ( 5||0'|| 2 + 5- 2 ||0|| 2 ), 

wmy\\ 2 <c(5n>\\ 2 +6~ 2 n\\ 2 ), [ ■ ' 

with C being constants not dependent on 5 and (/). 

Proof. We set a :— tt/2 as in the proof of Theorem 13.21 M is Hilbert-Schmidt 
since / = ftfr 1 = I + L + M + LM and L is Hilbert-Schmidt. 

We consider fi* at first. Following the proof of Theorem 13.21 L* can be 
written as 



oo 

L*f = J2Mx%J), (3.14) 

fe=0 



where <pk ■— 4>k — Xk an( l / G W. We show that L* is bounded on W 1 ' 2 {—a, a). 
We estimate the Hilbert-Schmidt norm of L* on W 1,2 (— a, a) with help of the 
orthonormal basis /„ :— Xn /Vl + n 2 . In fact, it suffices to estimate: 

oo oo 

Y J {(L*f n )',(L*f n )') = Y.YTn^^' j2 (3 - 15) 



where (recall (JHU) and (|37f|)) 

0n — ~~Tl£n Xn ^ti Xn S n Xn S n ^ti Xn ^n) 

+ ~ [Xn i 1 + ^ °n) - Z^Xn s n] ■ 



(3.16) 



Using the asymptotic properties of e n and the uniform boundedness of c„,s„ 
(see (|3.5p and (|3.7p . respectively) together with the normalization of x^, we 
conclude that \\<p' n \\ < C uniformly in n. Therefore (|3.15l) is finite. 

Using the same technique, we can show that the Hilbert-Schmidt norm of L* 
in W 2 ' 2 (—a,a) is finite. To this end we select the basis Xn + n 2 + n 4 , the 
rest is based on ||<^|| = 0{n) as n — > oo. 

Let us now establish the inequalities (|3.13l) . Consider <f> € W 1 ' 2 (—a,a), its 
basis decomposition <\> — Y^=o a nXn > an d the identity 

oo 

]T|na n | 2 = ||0'|| 2 . (3.17) 

n=0 
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Hence, 



oo 



ii(i»'ii 2 = E 



n 



(3.18) 



m,n— 



and having the explicit form of 4>' n , see (|3 . 16[) . we have to estimate several terms. 
We show the technique only for one term, the estimate of remaining terms is 
analogous. First, using the uniform boundedness of ||c„||, ||s„||, the asymptotics 
£„ = C(n _1 ) and the uniform boundedness of || II oo , it is easy to see that 



holds with some positive constant C . It remains to estimate the I -norm of a n 
by the I -norms of a n and na n (which equal ||0|| and \\4>'\\, respectively). This 
is rather algebraic: 



with any b, 5 € (0, 1). Here the first inequality follows by the generalized Holder 
inequality and the last one is a consequence of the Young inequality. The ex- 
ponent b is chosen in such a way that 2b > 1, so that the sum of n~ 2b (denoted 
by Cb) converges. If we put b = 2/3, we obtain the inequality in the claim. 

One can show, using the asymptotics (13.51) . that it follows from the normal- 
ization requirement ((t> n ,tp n ) — 1 that A n , the normalization constants of ip n , 
see (|2.16[) . satisfy A n = 1 + 0(n~ 1 ). Then the claims for O -1 and M can be 
derived in the same manner. 

To justify that and (fi^ 1 )* are bounded on W 1 ' 2 (—a,a) and W 2 ' 2 (— a, a), 
it suffices to realize that and fl* are invertible because they are invertible in 
L 2 (—a, a) and the inverse is bounded because of the form identity plus compact 
operator on considered Sobolev spaces. □ 

Corollary 3.7. Assume the hypotheses of Theorem ] 3. 61 Then h := iliJO -1 is 
a holomorphic family of operators of type (B) with respect to c± . The associated 
quadratic form th, in the sense of the representation theorem JUA Thm. VI. 2.1], 





<c b W\\ 2b U\\ 2(1 - b) 
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reads 



(V>(a) + (L*V)(a))(V>(a) + (M^)(o)) 



(^(-a) + (£*^)(-a)) ty(-a) + (M^)(-a)) 



(3.19) 



Dom(t h ) = W^(-a,a). 



Proof. The form t/, defined in (|3. 19|> is sectorial and closed due to the perturba- 
tion result [TH Thm. VI. 1.33], regarding u[tp] := th[t/j] — \\tp'\\ 2 as a perturbation 
of to[tjj} := ||^'|| 2 . Indeed, the inequalities (|3.13l) applied on u[tp] yield that u is 
io-bounded with to-bound 0. Therefore, due to the first representation theorem 
[T9l Thm. VI. 2.1], there is a unique m-sectorial operator associated with t/j. Let 
us denote it by h. Our objective is to show that h = h. 

Using the definition of h by the similarity transformation, i.e. h = fti/ft -1 , 
and the fact that H is associated to tn, we know that the domain of h are 
functions u such that, firstly, ft -1 ?/ € W 1,2 (— a, a) and, secondly, there exists 
w G L 2 {— a, a) such that 



(3.20) 



for all w such that ft*u G W 1,2 (—a, a). However, by Theorem l3~6l ft, ft*, ft -1 , 
(ft*) -1 are bounded on I^ 1,2 (— a, a) and it is easy to check that the identity 



i#(ft*v,ft 1 u)—th{v,u) 



(3.21) 



holds for all u, v € W ' (— a, a). Consequently, the operators h and /i indeed 
coincide. □ 

Remark 3.2. We remark that the boundedness of ft, ft*, ft -1 and (ft -1 )* in 
W 2 ' 2 (—a, a) was not used in the proof Corollarv l3.7l Nevertheless, this property 
is useful if we analyse the domain of h directly from the relation h = fti/ft -1 . It 
follows that Dom (h) consists of functions ip from W 2 ' 2 (—a, a) satisfying bound- 
ary conditions (ft~V)'(±a) + c± (ft - V)(±a) = 0. 



4 Closed formulae in PT-symmetric cases 

We present closed formulae of operators O, ft and h corresponding to H with 
special "PT-symmetric choice of boundary conditions, c± := ia, with a E R. 
This case has already been studied in a similar context in [2TJ [2U] , where the 
first formulae of the metric were given. We substantially generalize these 
results here. 

We essentially rely on the original idea of [3D] to "use the spectral theorem 
backward" to sum up the infinite series appearing in the definition of in (I2.2ip . 
The attempts to find ft as the square root of using the holomorphic and 
self-adjoint calculus are contained in [35| [34], however, only approximations 
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of the resulting similar self-adjoint operator h were found there. The main 
novelty of the present approach comes from the more general factorization (|2.22l) 
with (|3.1[) . which enables us to obtain exact results. Formulae contained in this 
section are obtained by tedious although straightforward calculations that we 
do not present entirely. 

Finally, we present the metric operator for H with general 'PT-symmetric 
boundary conditions, c± :— ia ± ft. In this case, the eigenvalues are no longer 
explicitly known, nevertheless, the experience from previous examples and for- 
mulation of partial differential equation together with a set of "boundary con- 
ditions" for the kernel of the integral operator provide the correct result. 



4.1 Reduction to finding a Neumann metric 

We start with the following fundamental result. 

Proposition 4.1. Let c± :~ ia, with a 6 R. Then the operator Q defined 
in (|2.21[) has the form 

e = j N + c'^e 1 + j N e 2 + J D e 3 , (4.i) 

where J L , with i € {D,N}, are defined in (|2.10[) . Co > 0, and 9i are integral 
operators with kernels: 



i(x,y) := ~e?( x y) sin (^(x - y)J 



h{x,y) := — [y-asgn(y-x)], (4.2) 
9 3 {x, y) := — {a - xy) - — x - — |1 - ia{y - x)\ sgn(y - x). 



is the metric operator for H , see Definition \2.3l if, and only if, a ^ k n for 
every n 6 N. 

Proof. Using the explicit form (I2.18[) of functions <fi n and the definition (|2.21[) 
of 9, we obtain 

oo 

e = E c nx%(x%, •> + cl (M<h, •) - xZixZ, •)) 

°° ^<2 00 00 ri2 

+ « 2 E w*Z<xZ' ■> + [a E irrt (x« . ->-i«E , •)• 

n=l 11 n=l n—1 

(4.3) 

Employing the operators J L andp,p* introduced in (|2.10[) and (|2.4[) . respectively, 



1G 



and relations (j2.5[) we obtain: 



n=0 

oo ^ 

+ «J iV p^^^,-) (4.4) 

n— 1 ft 

(OO OO \ 

n=l ™ n=l " / 

It follows from the functional calculus for self-adjoint operators that (j4.4|) can 
be written as 

6 = J N + Cl (Mh, •) - xZixS, •)) + a^pi-An)- 1 



By inserting the explicit integral kernels of the resolvents, see Section 12.21 we 
obtain the formula (f4~T|) with (|4.2|) . 

To ensure that such O represents as metric operator, we recall that the 
spectrum of H is always real, see Proposition 12 .41 Moreover, it is simple if, and 
only if, the condition in the last claim is satisfied. □ 

Remark 4.1. The formula ()4.1j) can be rewritten in terms of the operator J N 
only. Indeed, it is possible to show that 

J D =p*J N P (-A D )- 1 . (4.6) 

The final result is then 

e = j N + cle 1 + j N e 2 + P *j N e 4 , (4.7) 

where #4 := p{— Ar>)^ 1 9^ is an integral operator with kernel 



h(x, y) = — [y 2 {3 - iay) + 3x 2 (l - lay) + 2a 2 [l + ia(3x - y)] 
12a v 

- -a[ 2 - ia(y -x)\(y-x) sgn(y - x). 



(4.8) 



Note that the expression (|4.8|) is a result of a rather lengthy computation. 

Any metric operator for H in Proposition 14. II can be obtained by determin- 
ing J N for given constants C n . Thus we managed to transform the problem of 
constructing the metric operators for non-self-adjoint operator H to the prob- 
lem of constructing the metric operators J N for the Neumann Laplacian — Ajy- 
This significantly simplifies the problem, since — Ajv is self-adjoint and its met- 
ric operators are bounded, positive operators with bounded inverse commuting 
with — An. For instance, any bounded, uniformly positive function of — An 
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represents a metric operator. Moreover, it was shown in [33] that any J N can 
be approximated in the strong sense by a polynomial of I + A(-Ajv — A) , 
with A € p(— An). 

We consider two choices of constants C n in the following and we find final 
formulae for the corresponding metric operators. 

4.2 The constant-coefficients metric 

Let C 2 := 1 for every n > 0. Then J N = J D = I and the metric operator 
reads = I + K , where K is an integral operator with the kernel 

K(x,y) = - e ^ {x - y) sin (^(x-y)) + ^-(\y-x\- 2a) sgn(y - x) 

" V2 2a (4.9) 



2 



,2 







xy - a\y - x\) 



2a 

Formula (|4.9|) represents a remarkably elegant form for the metric operator 
found firstly in [21]l20]. 

4.3 The C operator 

Another choice of C n is motivated by the concept of C operator, see Defini- 
tion ED We want to find such that C 2 = I, where C = VQ. Since H is 
■P-self-adjoint, we have V<fi n = D n ip n with some numbers D n . Assuming the 
non-degeneracy condition a ^ k n for every n > 0, an explicit calculation shows 
that 

sin(2cra) k 2 — a 2 , 

Do = ^, D„ = (-!)» » 6 N. (4.10) 

The condition (TO) 2 = / then restricts C n from (|2~2Tj) to 

2 2|o|a 2 fc » c M /i ii\ 

In order to simplify the formulae, we consider only a £ (0, fci) in the following. 

Remark 4.2. As mentioned below (I2.21[) , any choice of C n can be interpreted 
as a sort of normalisation of <fi n . It is therefore interesting to notice that (|4. 1 1[) 
results into the symmetric normalization of <j) n and ip n when (</> n , ?/Vi) = 1 is 
required: 



M*) = \ -rr^ eiaae ~' iax > ^ = /, 2 " 2 (^^ ~ V^O*) 

V sm(2cra) ^/fc 2 - a 2 \ «n 

0o(*) = J^h;j aa e iax , <t>n{x) = kn (x%(x) + if t£{x) 

V sin(2aaj ^/fc 2 - a 2 \ fc n 

These expressions should be compared with the normalization of (|2.18[) (|2 . 1 9[) . 
standardly used in the present paper. The symmetric form of the "present 
normalization" indicates that the choice (|4. 1 1[) will lead to a simpler form of 
than (1431). 
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(4.12) 



Using (|4.1ip in the series (|2.10p . the operators J 1 can be determined by the 
functional calculus: 

jN = E in XnixZ, ■) + cl xZixZ, •> 

rZ OL 

71 — 11 

= {-A N ){-A N a 2 )- 1 + Cl x^(Xo W , •) 

= I + a\-A N -aY 1 +C 2 oX Z(x^-), 
00 , 2 
jD = y> fcn n, n x 

n=l n 

= (-A D )(-A D -a 2 )- 1 
= I + a 2 (-A D -a 2 )-\ 

A direct (but very tedious) way how to derive the metric O for the choice (|4. 1 1[) 
is to express the resolvents of the Dirichlet and Neumann Laplacians from the 
ultimate expressions in (|4.12l) by means of the Green's functions (|2.6p and com- 
pose them with the operators 9i in (|4.1[) . 

However, a more clever way how to proceed is to come back to the operator 
form (|4.5[) and perform first some algebraic manipulations with the interme- 
diate expressions appearing in (|4.12l) . First, we clearly have J D (— A^) _1 = 
(—A D - a 2 ) -1 . Second, employing (j2l| and the identity (-A A r)(-A^)~ 1 = 
I-Xo (Xor), we check 



2 > 



[jV(-Ai)- 1 ]* - P (-A D - a 2 )~\ [jKpi-An)- 1 ] * = p*(-A N - a 

Finally, again using (|2.4j) . we verify the intertwining relation [p (-Afl-a 2 ) -1 ]* = 
p*{— Ajy — a 2 ) -1 . Summing up, with our choice (|4.11[) . formula (|4.5p simplifies 
to 

6 = / + Cg M<t>o, ■) + a 2 (-A N - a 2 )- 1 + a 2 (-A D - a 2 )' 1 
+ ap(-A D -a 2 )- 1 +ap*(-A N -a 2 )- 1 . 



Now it is easy to substitute (|2 .6[) and after elementary manipulations to conclude 
with = 1 + K, where K is an integral operator with the kernel 

JC(x,y) =ae- ia{y -' x) [tan(cm) - isgn(y - x)]. (4.14) 

The operator C can be found easily by composing V and O. We finally arrive 
at the formula C = V + L, where L is an integral operator with the kernel 

C(x,y) = ae-' ia{v+x) [tan(aa) - isgn(y + x)] . (4.15) 
4.4 The similar self-adjoint operator 

Next we present an example of operator J7, defined in (13.11) with e„ := Xn 1 that 
will be used to find the similar self- adjoint operator h from (|1.3j) . We recall 
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that the similarity transformation is invertible if all the eigenvalues of H are 
simple, which is ensured by the condition a ^ k n for every n G N. We will 
actually search for the quadratic form associated to h for which we have the 
result in Corollary 13.71 

We follow the analogous strategy to obtain formula for f2 as in the proof of 
Proposition 14. II The definition of f2 with e n := Xn leads to the sum: 



n = •> + £ xM, ±xZ(xZ, ■ 

i i "Vi 

n— 1 71=1 
oo 



n=l n 

I + X${4>0, •) - xZixZ, •) + ap(-A D )-\ 



(4.16) 



where we have used identities (|2.5[) . In the same manner, we obtain the result 
for the inverse fl^ 1 : 

^ fc « - a ,tt fc « - a (4.17) 

= I + V>o<X^, •) + a 2 (-A N - a 2 )- 1 - ap*(-A N - a 2 )' 1 . 

The operators L, M appearing in the expressions for fl = I+L and O -1 = I+M 
are, as expected, integral operators with the kernels C, M. that can be easily 
obtained using formulae for the Neumann and Dirichlet resolvents (|2.6|) - (|2.7[) : 



C(x, y) = ¥L[y- a sgn(y - x)] + (e~ ia ^ - 1 

np ia(a— x) n 

M(x, y) = . j- . - - e- ia ^[cot(2aa) ~ isgn(y - x)] (4.18) 



ae 



-ia(x+y) 



2 sin(2aa) 

To find the similar self-adjoint operator (jl.3l) . we start from the quadratic 
form (|3.19l) . Inserting (|4.18l) into the latter and performing several integrations 
by parts with noticing that LM = — L — M and (Mtp)' = —laM tp — iaip results 
in: 

t h [i>] = \\ip'\\ 2 + a?\(xS,Tp)\ 2 - (4-19) 
The corresponding operator h reads: 

Dom(ft) = {ip e W 2 - 2 (-a,a) : ip'(±a) = 0} . 

We remark that h is a rank one perturbation of the Neumann Laplacian. The 
eigenfunctions of h are Xn with Xo corresponding to the eigenvalue a 2 . 

It is interesting to compare the spectra of H and h for a = k n , i.e. in the 
points where the spectra are degenerate and similarity transformation breaks 
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down because the operator fl is not invertible. fc^ is an eigenvalue with the 
algebraic multiplicity two for both H and h. However, the geometric multiplicity 
differs: it is one for H and two for h. 

The form of h also explains the origin of the peculiar a-dependence of the 
eigenvalues of H (which are all constant except for \o(a) = a 2 ). In fact, it is 
the nature of the rank one perturbation to leave all the Neumann eigenvalues 
untouched except for the lowest one that is driven to the a 2 behavior. 

4.5 More general boundary conditions 

Finally, we consider the general 'PT-symmetric boundary conditions c± :— ia± 
f3, with a !( 8 G R. We start with formal considerations. The 9-self-adjointness 
of H can be expressed in the following way. We take the advantage of the 
realization of O = I + K, which we insert into QHip — H*Qip, ip £ Dom (H). A 
formal interchange of differentiation with integration and integration by parts 
yield following problem that we can understand in distributional sense: 

(%-tf)K(x,y)=0, (4.21) 
d y K(x, ±a) + (ia ± 0)K{x, ±a) = 0. (4.22) 

Moreover, Qip must belong to Dom (H*), from which we have a condition 

d x K(±a, y) + (-ia ± /3)/C(±a, y) = 2ia% T a). (4.23) 

Here S denotes the Dirac delta function. 

Already presented examples of for f3 = satisfy these requirements, par- 
ticularly K. solves the wave equation (I4.21j) . The kernel (|4.14p . corresponding 
to the simpler form of presented metric operators, is a function of x — y only. 
Inspired by this, we find the solution of the wave equation 

JC(x,y) = e ia (*-v)-f>\*-v\ [c + ia sgn(x - y)] , c e R, (4.24) 

that satisfies the "boundary conditions" (|4.22p and (|4.23|) as well. The one 
parametric family of solutions (I4.24j) of (|4.21l) - (|4.23l) demonstrates the known 
non-uniqueness of solutions to this problem. We also remark that c can be taken 
as a or a dependent as well. 

The positivity of O is ensured if the norm of K is smaller than 1. This can 
be estimated by the Hilbert-Schmidt norm of K which is explicitly computable: 

l|if|| H s = (c 2 +« 2 ) P 2/g2 • (4.25) 

Consequently, the positivity of & can be achieved by several ways, e.g., if a is 
small; or if /3 is positive and large; or \c\ and \a\ are small. In any of the regimes, 
the formal manipulations above are justified. 



21 



5 Bounded perturbations 



In this section we show that results of Section [3] remain valid if we consider a 
bounded perturbation V of H. 

Firstly we remark that the perturbation result [TTJ Thm. XIX 2.7] guar- 
antees that H + V remains a discrete spectral operator. That is, if all the 
eigenvalues of H + V are simple, then the metric operator O exists. We show 
that the claim of Theorem l3.2l is valid for H + V as well. The rest of the results 
from Section [3] then follows straightforwardly. 

Our approach is to use analytic perturbation theory for the operator h := 
fllin^ 1 that is perturbed by a bounded operator f2VT2 _1 . We denote by £„, r/ n 
the eigenfunctions of H + V and H* + V* , respectively. Let e„ be elements of 
any orthonormal basis in %. 

Theorem 5.1. Let all the eigenvalues of H be simple and let V be a bounded 
operator. If all eigenvalues of H + V are simple, then fiy = J2^=o e n(Vn, ■), 
i.e. fly '■ £n I - * e„, can be expressed as 

n v = U + L, (5.1) 

where U is a unitary operator and L is a Hilbert- Schmidt operator. 

Proof. As in the proof of Theorem 13.21 without loss of generality, we restrict 
ourselves to e n := Xn an d we show that fiy = I + L with L being Hilbert- 
Schmidt. We consider the normal operator h := flHfl^ 1 and we perturb it by 
v := ilVfl^ 1 . More specifically, we construct h{e) := h + ev forming a holomor- 
phic family of type (A) with respect to the parameter e. We denote by /Lt n (e), 
/Un(e) the eigenvalues and by £n(ff), fjn(s) the corresponding eigenfunctions of 
h(e) and of h(e)* respectively. h(0), h(0)* are normal, therefore the eigenfunc- 
tions £ rl (0) and rj n (0) form orthonormal bases. In fact, with our choice of e„, 
£ n (0) = r?„(0) = x». 

We construct operator ft : £, n (l) i— > Xn an( i we wm show that Cl = I + L, 
where L is Hilbert-Schmidt. Qy is the composition of fi and f2 and the claim 
then follows easily using of the fact that Hilbert-Schmidt operators are a *-both- 
sided ideal. 

The distance of /i„(0) and jtt n (l) can be at most \\v\\. Since we know the 
asymptotics of fx n (0) = A rl , see (|3.5j) . it is clear that there exists A^o such that 
for all n > No, |/i„+i(l) — fj, n (l)\ > n holds. Moreover, for such n the radius of 
convergence of perturbation series for eigenvalues and eigenfunctions is larger 
than 1. Thus, we have 

oo 
3=1 

(i) 

We estimate the norms of fjn using the analytic perturbation theory: 

7^11 < ^- i \\(h(oy-E)- l (v*(h(oy-E)-y x %\\dE 



3 r J 



1 / V+ 1 v 



(5.3) 
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where T n is a circle around /x n (0) of radius n/2 and the constant c does not 
depend on n. We define Ni as such that Ni > No and c/iVi < f . 

We prove that Cl has the desired form by showing that the adjoint Cl* = 
Sr=o Vn(^)(Xn > •) can be written as ft* = I + L 



JVi-1 



oo oo 



L^, where 



T* 



n =iVl j = l 



(5.4) 



and L^jy and are Hilbert-Schmidt. The decomposition of fl* follows im- 
mediately if we consider the expansions (|5.2p for n > iVi and rewrite fj n (l) = 
Xr^ + (^n(l) — Xn) f° r n < N\. Ljy is a finite rank operator therefore it is 
obviously Hilbert-Schmidt. is bounded and the defining sum is absolutely 
convergent since 



oo oo 



E Eii^iikx-^m E E© E 

n=Wi j'=2 n=iVi 3=2 n=iVi 



nc 



E II^IIKx^M^ 

n— A^i 



oo _^ 
\ E 



E i<xMi 2 < C |Hi 

n=Ni 



\ 



OO ^ 

T1 2 



(5.5) 



n=N\ 



Finally we estimate the Hilbert-Schmidt norm of L* 



oo oo 



oo oo 



^ ^ ( ^ ^ ^ ^ Vrii iXmi Xp ) i ^ ^ ^ ^ ^ (Xn i Xp ) 

3 oo / x 2 



oo oo 



^EE 

p=N x i=l 



E 



^E 



p — c 



< oo. 



This concludes the proof of the theorem. 



(5.6) 



□ 



Remark 5.1 (General Sturm-Liouville operators). Let us conclude this section 
by a remark on how to extend the previous result on bounded perturbations V 
for the operator H in the general form 



on 



Htp := -(pip')' +Vijj 
subject to the boundary conditions 

p(±a)%p' (±a) + c±ip(±a) 



L 2 (-a,a) , 



0. 



(5.7) 



Assuming merely that p is a bounded and uniformly positive function, i.e., 
there exists a positive constant C such that C _1 < p(x) < C for all x € (—a, a), 
the operator can be defined (cf. |8J Corol. 4.4.3]) as an m-sectorial operator 
associated with a closed sectorial form with domain W 1,2 (— a, a). If, in addition, 
we assume that p G W 1,00 (—a,a), then it is possible to check that the domain 
of H consists of functions ip from the Sobolev space W 2,2 (— a, a) satisfying (|5.7p . 
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Now, let us strengthen the regularity hypothesis to p G W 2 -°°(—a, a) and in- 
troduce the unitary (Liouville) transformation U : L 2 (—a,a) — > L 2 (f(—a),f(a)) 

by 

U- l <t>:= /9" 1/4 0o/, where fix) := [ — %= . 

Then it is straightforward to check that the unitarily equivalent operator H := 
UHU- 1 on L 2 {f(~a),f(a)) satisfies 

H(j) = -<f>" +V<f>+ W<f>, 
Dom(H) = {<(> eW 2 > 2 (f(-a),f(a)) : 0'(±/(o)) + £ t ^(±/(a)) = 0}, 

where V :=UVU- Y and 

^(ia) 1 ^ 4/>(±o)V2' ■ 16 p 7 7 ' 

In this way, we have transformed the second-order perturbation represented by p 
into a bounded potential 14 7 and modified boundary conditions. Theorem 15.11 
applies to H and, as a consequence of the unitary transform U, to H as well. 



6 Conclusions 

In this article, we investigated the structure of similarity transformations f2 and 
metric operators for Sturm-Liouville operators with separated, Robin-type 
boundary conditions. The main result is that 51 and can be expressed as a 
sum of the identity and an integral Hilbert-Schmidt operator. 

We would like to emphasize that this not always the case for other types 
of operators, see, e.g., [2 [32J, [24j [14], where is a sum of the identity and 
a bounded non-compact operator. The latter is a composition of the parity 
and the multiplication by sign function. Moreover, corresponding similarity 
transformations map (non-self-adjoint) point interactions to (self-adjoint) point 
interactions, which is not typically the case for operators studied here. This 
is illustrated in the example of ■pT-symmetric boundary conditions where the 
equivalent self-adjoint operator is not a point interaction but rather a rank one 
perturbation of the Neumann Laplacian. 

In this work we considered the separated boundary conditions only. Nonethe- 
less, the analogous results are expected to be valid for all strongly regular bound- 
ary conditions. 

As the proofs of the results show, the crucial property is the asymptotics of 
eigenvalues, i.e. separation distance of eigenvalues tends to infinity, that is used 
for the proof of the existence of similarity transformations [11] . Recent results on 
basis properties for perturbations of harmonic oscillator type operators [TJ [30l 13] 
give a possibility to investigate the structure of similarity transformation in these 
cases as well. Another step is to extend the results e.g. on Hill operators, where 
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a criterion on being spectral operator of scalar type has been obtained in [IB") 
and recently extended in [ID) . 

On the other hand, the structure of similarity transformations for operators 
with continuous spectrum as well as for multidimensional Schrodinger operators 
is almost unexplored and constitutes thus a challenging open problem. 

We illustrated the results by an example of 'PT-symmetric boundary condi- 
tions, where we found all the studied objects in a closed formula form, which is 
hardly the case in more general situations. However, in general, we may search 
for approximations of O or 9, typically applying the analytic perturbation the- 
ory to find perturbation series for eigenvalues and eigenfunctions of H to certain 
order k. For instance, we perturb the parameters c± in boundary conditions 
by small e. As a result we find an approximation h app of the similar operator 
h with resolvents satisfying \\(h — zy 1 — (h app — z)^ 1 ]] < Ce k . An extensive 
discussion and example of such construction can be found in [34] . The same 
remark is appropriate for small perturbations by bounded operator discussed in 
Section [5] 
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